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Abstract 

The large-x behavior of the transverse-momentum dependent quark distributions is analyzed in 
the factorization-inspired perturbative QCD framework, particularly for the naive time-reversal- 
odd quark Sivers function which is responsible for the single transverse-spin asymmetries in various 
semi-inclusive hard processes. By examining the dominant hard gluon exchange Feynman diagrams, 
and using the resulting power counting rule, we find that the Sivers function has power behavior 
(1 — x) 4 at x — > 1, which is one power of (1 — x) suppressed relative to the unpolarized quark 
distribution. These power-counting results provide important guidelines for the parameterization 
of quark distributions and quark-gluon correlations. 
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I. INTRODUCTION 



Single transverse-spin asymmetries (SSA) have a long history, starting from the observa- 
tion of large SSAs in hadron production in nucleon-nucleon scattering in late 70s and 80s 
PJ. Since initial or final state phases are required to produce these T — odd observables, 
SSAs provide a unique window into quantum chromodynamics (QCD) at the amplitude level 
as well as the role of quark orbital angular momentum in the wave functions of hadrons. 
Recent experimental observations of sizeable SSAs in hard scattering reactions such as single- 
inclusive deep inelastic scattering £p — > £'tt + X have greatly motivated new studies of the 
underlying mechanisms in QCD. These experiments include semi-inclusive deep inelastic 
scattering at HERMES at DESY @, COMPASS at CERN 0, and Jlab 4], and hadron 
production in nucleon-nucleon scattering at RHIC 0,0,0]. Large SSAs have been observed 
in semi-inclusive production of hadrons in DIS for transversely polarized proton target at 
HERMES 0, and single inclusive hadron production in the forward direction in polarized 
proton-proton scattering at RHIC 0,0]. 

On the theory side, two mechanisms have been proposed in the QCD framework to explain 
these large SSAs in hard scattering processes. One is based on the QCD collinear factor- 
ization where the asymmetries arise from the higher-twist quark-gluon correlation effects 
(Efremov-Teryev-Qiu-Sterman mechanism) 0, 0. Another approach explicitly takes into 
account the effects coming from the intrinsic transverse momentum of partons in hadrons. 
For example, the Sivers function was proposed in |l0| to explain the SSA phenomena in 
hadronic reactions, where intrinsic transverse momentum plays an important role. 

In the last few years, there has been an intensive theoretical development of transverse mo- 
mentum dependent (TMD) parton distributions and their roles in semi-inclusive processes 
such as semi-inclusive deep inelastic scattering (SIDIS) and the small transverse momen- 
tum Drell-Yan process. The gauge-invariant properties [ill E2I. Il3l llil of the TMD parton 
distributions and the relevant factorization formalism |l5l Il6l Il7l Il8l| have been studied 
thoroughly. For example, the Sivers effect in SIDIS has been shown to arise from the in- 
terference of amplitudes differing by one unit of quark orbital angular momentum and the 



fact that these amplitudes have different final-state phases |ll[. The phases arise from the 
Wilson-line associated with the struck quark as required by gauge invariance 12||. The SSA 
reverses sign in Drell-Yan reactions because the phases in the Drell-Yan reaction arises from 
initial-state rather than final-state interactions 11, Remarkably, the SSA effect in these 
reactions is leading twist; i.e., it survives in the Bj or ken- scaling limit. Moreover, it was 
recently shown that the above two mechanisms for SSAs are unified for physical processes 
in the kinematical region where both apply [l9| . 

There has also been a number of phenomenological studies of the experimental data. 
Model-dependent parameterizations of the relevant non-perturbative parton distributions 
(twist-3 quark- gluo n correlation or the TMD quark distributions) have been adopted to fit 



to the data |9|, |20j, |21|, |22j, 23] . In these studies it has been implicitly assumed the Sivers 
function is suppressed at large x relative to the unpolarized quark distributions (21, 22||. In 
this paper, we will provide an argument for this suppression based on power-counting of the 
leading diagrams in perturbative QCD. We will utilize the generalized power counting rule 
and adopt a perturbative analysis of the structure function at large x. 

The large-x behavior of both the polarized and unpolarized parton distributions have 
been studied ji^, H^, [2^, 27, 2^] in PQCD. A generic factorization has recently been used 



to justify the power counting rule by relating parton distributions at large-x to the quark 
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distribution amplitudes of hadrons 29]. So far, the power counting results have been worked 
out for the unpolarized and longitudinal polarized quark distributions. In the present study, 
we will extend this analysis to other leading-order TMD quark distributions, including the 
naive time-reversal-odd quark Sivers function which is responsible for the SSAs in various 
semi-inclusive hard processes. 

It is important to note that the x — > 1 regime where the struck quark has nearly all of 
the light-cone momentum of its parent hadron involves dynamics far-off the mass shell: the 

Feynman virtuality of the struck quark becomes highly space-like: kj^—m 2 ~ — - , where 
k± and A4 are the transverse momentum and invariant mass of the spectator system. Thus 
we can use perturbative QCD to analyze the large-x behavior of parton distributions since 
the internal propagators in the relevant Feynman diagrams scale as 1/(1 — x). This behavior 
leads to a power counting rule. This is because more partons in hadron's wave function means 
more propagators in the scattering amplitudes, and more suppression for the contribution to 
the parton distributions. Thus the parton distributions at large-x depend on the number of 
partons in the Fock state wave function of the hadron. In particular, the valence Fock state 
with the minimum number of constituents will dominate the quark distribution function 
at large-x. For example, the proton structure function will be dominated by its three- 
quark Fock states, which can be further classified according to its quark orbital angular 
momentum projection: L z = 0, \L Z \ = 1, or \L Z \ = 2 30]. Since nonzero quark orbital 
angular momentum light-cone wave function normally introduces additional suppression of 
(1 — x), we will consider in this paper only L z = and \L Z \ = 1 Fock states contributions. 
The \L Z \ = 1 state is needed because some of the TMD quark distributions involve the 
interference between L z = and L z — 1 states (HI HH Hl| (see also the discussions below). 

As is the case of the nucleon form factors (Dirac and Pauli form factors) j^, the trans- 
verse momentum dependent quark distributions can be calculated from the overlap of the 
light-cone wave functions of three-quark Fock states jTH HO, Hl| . As we shall demonstrate, 
the large-x power counting for the TMD parton distributions can be obtained in a similar 
manner. For example, we know that the unpolarized quark distribution has power counting 
of (1 — x) 3 at large-x |24j, which comes from the quark orbital angular momentum projection 
L z = Fock states contribution, whereas the contribution from the overlap of two light-cone 
functions for \L Z \ = 1 states is suppressed by (1 — x) 2 j2^, |2(| 27, 28||. On the other hand, 
since the Sivers function depends on the interference between L z = and \L Z \ = 1 states, 
simple counting suggests that that the Sivers function will have the leading power of (1 — x) 4 . 
The detailed calculations in this paper support this intuitive argument. 

The remainder of the paper is organized as follows. In Sec. II, we present our analysis of 
the leading-order TMD quark distributions at large-x, where we discuss the power counting 
results for the k±-eveia, k±-odd, and naive time-reversal-odd quark distributions respectively. 
We will also derive the power counting results for the integrated parton distributions at 
leading-twist and sub- leading- twist. We summarize our results in Sec. III. 



II. TRANSVERSE-MOMENTUM DEPENDENT QUARK DISTRIBUTIONS AT 
LARGE-X 



The TMD quark distributions can be defined through the following matrix: 

di~d 2 U 



M 



ai3 



(27T) 



E?-P++i£±-fc± /pg 



PS 



(1) 
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FIG. 1: Typical Feynman diagram contributing to the large-x quark distribution in nucleon. The 
blobs at the left and right sides represent the three-quark light-cone wave function distribution 
amplitudes of the nucleon. 

where the vector P = (P + , 0~, 0j_) is along the momentum direction of the proton, 5* is the 
polarization vector, and & v (£) is defined as 

^(0 = ^(005^(0 , (2) 

with the gauge link C v (oo; = ex P ( — W d^v ■ A(\v + 0)- In this paper, we will study 
the TMD quark distributions for the semi-inclusive DIS processes, thus the above gauge 
link goes to +oo. Our results can be simply extended to the TMD quark distributions for 
the Drell-Yan process, where the naive time-reversal-odd TMDs will have an opposite sign 
(Tl . 12l . 13| . In the following analysis, no light-cone singularity will be present, thus we can 



choose the vector v to be a light-cone vector n = (0 + , n~, 0j_) with n ■ P = 1. 

The leading order expansion of the matrix M. contains eight quark distributions |3c 
Among the eight distributions, three are the so-called fcj_-even TMD quark distributions: 
q(x, k±), Aq L (x, k±), and 5qr(x, k±), which correspond to the unpolarized, longitudinal po- 
larized, and transversity distributions, respectively. These distributions will lead to the 
three leading-twist integrated quark distributions |34| after integrating over transverse mo- 
mentum. The other five distributions are k±-odd, and vanish when k± are integrated. Two 
of them, qr(x, kj ) (Sivers) and 5q(x, k±) (Boer-Mulders), are odd under naive time-reversal 
transformatio n 13511 . The notations for these distributions follow Ref. [30], which are different 
than those in |33|. However, their definitions are identical. 

In this paper we are interested in studying the large-x behavior of these TMD quark 
distributions. We will disregard the k± dependence, and further choose k± ^> Aqcd in order 
to avoid the infrared divergence associated with low transverse momentum limit. A typical 
Feynman diagram contributing to large-x quark distributions is shown in Fig. 1. At this 
order, we can write down an inspired factorization formula for the the parton distributions 
in terms of the distribution amplitudes of the nucleon (jR |26| , 

pi 7 \ f d?ki±cPk2± dz\dz2 , , , , ~, f r , ir , 

f(x,k±) = J ^ 2 7i) 6 ^—o{k ± + k 1± + k2±)d{z 1 + z 2 + x-l) J [dyiWdy^ 

Qiy^'iy'jHiyuy^h^Zi) , (3) 

where the outside integral represents the phase space integrals for the final state two quarks 
going through the cut line, with momenta: kj = (ziP + , k^ , k{±) (i = 1,2). The inside 
integral measure [dyi\ is defined as [dyi\ = dy\dy2dy$5{l — 2/1 — 2/2 — 2/3)5 and the yi are the 
momentum fractions of the proton carried by the quarks in the light cone wave functions, 
i.e., Pi = y^P and •p\ = y\P in Fig. 1. Here / represents any of the leading order TMD quark 
distributions. $ and $' represent the quark distribution amplitudes of nucleon at the left 
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and right sides of the cut line, respectively. They can be the leading-twist or higher-twist 
distribution amplitudes, depending on the quark orbital angular momentum projection along 
^-direction. We list the distribution amplitudes we will use in this paper in the Appendix 
for reference. 7i represents the hard part which can be calculated from the perturbative 
Feynman diagram like Fig. 1. 

We notice that, in the above equation the phase space integral for ki and k 2 are strongly 
constrained in the limit of x — > 1, because of momentum conservation, z\ + z 2 = 1 — x. 
We can factor out the (1 — x) dependence of this phase space integral, taking the following 
parameterizations: z\ — <x(l — x), z 2 = (3(1 — x), 

f dz 1 dz 2 1N 1 f dad[3 n , . A . 

—^5(z 1 + z 2 + x-l) = - —f6(l-a-0). 4 

J z x z 2 l-x J a/3 

This leads to an overall enhancement of 1/(1 — x). After factoring this out, the remaining 
measure of the phase space integral (dad[3) does not contain any additional factors of (1—x). 

Additional (1 — x) factors can come from the hard amplitude TC, but these depend on the 
structure of the relevant tree diagrams. Since the hard propagators each contain a (1 — x) 
factor, the least number of active particles involved in the hard process lead to the least 
suppression. Thus the leading contribution to the quark distributions at large-x is dominated 
by the leading component in the hadron's Fock state expansion. For the nucleon, the three- 
quark Fock state components will dominate the quark distributions, while for pion it will be 
the quark-antiquark pair states. In the following we will study the large-x power counting 
for the above mentioned TMD quark distributions, including the three fc^-even ones: q, 
Agz,, 5qr, and four k±-odd ones: 5qL, Ag^, <?t, and 5q. For 5q' T , its analysis will involve 
much more complicated diagrams, and we will not discussed this in the present paper. 



A. k±-even quark distributions 

The unpolarized quark distribution is defined as 

k±) = \ J ^r 6 ^^^ ( P 1^(07^(0) | P) ■ (5) 

The large-x power counting for this distribution function has been studied in the literature 
jiil lij^ . In the following, we will repeat these arguments as guideline for the analysis of 
other quark distributions. 

We calculate the above matrix element in the proton helicity basis, 

q{x > h±) = \j ^WF e ' ixrP++ ^ ±k± ^ (< P5 *l a l p5 *> + ( ps *i\ d \ ps *i)) > ( 6 ) 

where the operator O is defined as O = !P' t ,(^)7 + !P r 1 ,(0). This operator is chiral-even, and 
conserves the quark helicity in the partonic scattering matrix elements, and so that the 
dominant contributions come from the leading Fock state wave function (L z = 0) at both 
sides of the cut in Fig. 1. The matrix element (PS Z ^\(D\PS Z ^) will have the contributions 
from the following quark spin configurations: |J,f and where in the first one the probing 
quark has the same helicity as the proton and in the second case it is opposite. If the probing 
quark's spin is parallel to the proton spin, the two spectator quarks will form a scalar. 
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(a) (b) (c) 



i -^-J 



(rf) (e) 

FIG. 2: Leading diagrams contributing to the k±-even quark distributions at large-x: left half of 
the relevant diagrams are shown. The contributions will the amplitudes square of these diagrams, 
including the interference between them. These diagrams also contribute to the k±-odd and naive 
time-reversal-odd TMD quark distributions. 



In the following, we are only interested in obtaining the power counting for the quark 
distributions, and the explicit dependence on the distribution amplitudes will not be dis- 
cussed. According to the factorization formula Eq. Q and the reduced integral Eq. (jlj), the 
unpolarized quark distribution depends on the leading-twist distribution amplitudes of the 
nucleon, 

gfofcOUioc— !— / ^^Mvi^MHiviiVi^'PAi-x)) ■ ( 7 ) 

1 — x J ap 

The power counting of the hard factor 7i can be evaluated from the partonic scattering 
matrix elements, which include a set of propagators and traces of the Dirac matrices. As 
mentioned above, the propagators are far off-shell in the limit of x — > 1, which will lead to 
suppression in terms of (1 — x). For example, One of the gluon propagator in Fig. 1 goes 
like, 

1 _ 1 1 1 -x 

at large x. In the above expression, we have omitted all higher order terms suppressed by 
(1 — x). (k±) represents a typical momentum scale in order of transverse momentum k±. 
Besides the propagators, the traces of the Dirac matrix contains (1 — x) factors as well, 
which will depend on the spin structure of the quarks in the scattering matrix elements. 
For example, the traces of the Dirac matrices for the scattering Tit - *TIT i n Fig- 1 wm 
contribute to the matrix element (PS z ^\0\PS z ^) as 

in the leading power, where the probing quark has the same helicity as the proton. If the 
probing quark has an opposite helicity as the proton, e.g., in the spin structure the 
Dirac trace for the diagram of Fig. 1 vanishes. We have checked all other diagrams, and 
found that the Dirac traces for those diagrams with spin-one configuration for the spectator 
quarks (ft or II f° r the two spectator quarks) either vanish or are suppressed by at least 
(1 — x) 2 as compared to the scalar configuration (|| or []). This property has long been 
noticed in the literature 0, HE H3|- m Fig. 2, we showed all leading diagrams for the 
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spin structure of fit, where only left half sides of the relevant diagrams are shown. The 
contributions will be the amplitudes square of these diagrams, including their interferences. 

In summary, the leading contribution to the matrix element in Eq. © comes from the 
quark spin structure with probing quark's helicity equal to the nucleon's helicity. If they 
differ, the contribution will be suppressed by (1 — x) 2 . This is to say, the quark distribution 
is dominated by the quark spin parallel to the nucleon spin, and the quark spin anti-parallel 
distribution will be suppressed by (1 — x) 2 . 

The final power counting result will depend on the above factors, including the Dirac 
matrix traces, the power counting of the propagators, and the phase spaces integrals. For 
example, the diagram in Fig. 1 contains eight propagators, with overall power counting, 

~ y 3 (l - yi) 2 y' 3 (l - y[) 2 

Combining the above with the contribution from the Dirac matrix traces in Eq. © and the 
phase spaces integral factor in Eq. (j3J), we find that the contribution of this diagram to the 
unpolarized quark distribution is 

g(x,fc ± )|^ 1 ~(l-x) 3 . (ll) 

All other diagrams in Fig. 2 will contribute the same power to the unpolarized quark distri- 
bution. For example, the amplitude square of Fig. 2(d) contribute a power of (1 — x)~ 2 from 
the Dirac matrix traces, a power of (1 — x) 6 from the propagators, plus a power of (1 — x)~ l 
from the phase space integral Eq. (@J), which leads to a power of (1 — x) 3 contribution to the 
quark distribution. 

The longitudinal polarized quark distribution can be analyzed accordingly, which is de- 
fined through the following matrix element, 

Aq L (x, k 1 _) = \J ^^e-^-P^-k, ( PSz |^(O 7 + 75^(0)| PS,) . (12) 

Again, if we calculate in the proton helicity states, we will obtain 

Aq L (x,k ± ) = ^J ^^ e -^- p++ *^ , (13) 

where the operator Ol is defined as Ol = ^(£)7 + 75^(0)- We can interpret longitudinal 
polarized quark distribution as the quark spin parallel to the nucleon spin distribution 
minus the antiparallel distribution. According to the above analysis for the unpolarized 
quark distribution, we know that the quark spin parallel to the nucleon spin distribution 
dominates over the antiparallel distribution, and the latter is suppressed by an extra factor 
of (1 — x) 2 . In conclusion, we will obtain the same power behavior for the longitudinal 
polarized quark distribution as the unpolarized quark distribution, 

Ag L (x,fc ± )|^! ~ (1-x) 3 . (14) 

The transversity distribution for the quarks can be analyzed in the same way. It is defined 

as 
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(10) 



Here, the proton is transversely polarized. We can choose the polarization vector along 
x-direction, and the polarization states can be constructed from the proton helicity states, 



|P5 st > = \ (\PS z1 ) + \PS zl }) , \PS xl ) = 4= (\PS z1 ) - \PS zl }) 



V2 



V2 



(16) 



Substituting the above into the definition of the transversity distribution, we will obtain 



d£ d 2 £± -ix£-P++i£_ L -k ± 1 

(2tt) 3 



-l(P3^\O t \PS 



PS 



zl 



(17) 

where the operator O t is defined as O t = & v (£ , 0, &±)7 + 7" L 75^ r 1) (0). From this equation, we 
see that the quark transversity distribution depends on the matrix elements with hadron 
helicity flip. On the other hand, because the operator Ot is chiral-odd, it changes the 
quark helicity in the partonic scattering process as well. If we keep the leading Fock state 
contribution, the matrix element in the bracket of the above equation will reduce to 



PS z] \O t \PS 



A 



PS 



A 



o t 



PS 



A 



(18) 



where the subscripts ±^ represents the total quark helicity in the three-quark wave function 
used in the calculations. From this equation, we can easily see that it will be the same set 
of diagrams in Fig. 2 contributing to the transversity quark distributions. The same power 
counting results will be obtained, 



8q T (x, k 



(1-x) 



(19) 



From the above analysis, all the three k±-even quark distributions have the same power 
behavior at large-x, which is certainly consistent with the inequality condition for them [36] . 



B. fcx-odd and naive time- reversal- even quark distributions 

In this subsection, we will study two fcj_-odd but naive time-reversal-even TMD quark 
distributions: AqT and 5qT, which represent the longitudinal polarized quark distribution 
in a transversely polarized proton and the transversely polarized quark distribution in a 
longitudinal polarized proton, respectively. AqT can be calculated from the following matrix 
element, 

where the operator Ol as defined above. Following the above calculation for the transversity 
distribution, we choose the transverse polarization vector along the x-direction, and the 
above equation can be reduced to 



Or 



PS 



(21) 

Ol is a chiral-even operator, and it conserves the quark helicity. On the other hand, the 
above matrix element has hadron helicity flip, thus the total quark helicity and the hadron 



8 



helicity will mismatch on either side of the above matrix element. If the total quark helicity 
and the proton helicity is mismatching, the wave function for the three-quark state must 
have nonzero quark orbital angular momentum. It is the interference between the L z = 
and \L Z \ = 1 states contributing to the TMD quark distribution Aq T . 

In order to proceed, we further decompose the proton spin state into the Fock states 
containing L z = and \L Z \ = 1. For example, 

\PS z1 ) = \PS z1 ) 1/2 +\PS z1 )„ 1/2 , (22) 

where the subscript ±1/2 denotes the total quark helicity. The first term in the above 
equation represents the L z = state, while second one for the L z = 1 state. The wave 
function parameterizations for these states have been given in Eq. (jAlj) . Similarly, for 

\PS z i) we have, 

\PS zl ) = \PS zj )- 1/2 + \PS zl ) 1/2 , (23) 

where the first one is for L z = 0, and the second one for L z = —1. Because the partonic 
matrix element conserves the quark helicity, in the calculation of the matrix element of 
Eq. (|21|). the quarks helicities will remain the same at the left and right sides of the cut 
line in the Feynman diagram like Fig. 1. From the experience in the last subsection, we 
know that the partonic processes where the two spectator quarks have opposite helicities 
dominate the quark distributions at large-x. So, for the leading contributions we will have 
two typical partonic processes: TIT - ►TIT with total quark helicity 1/2 and |T! — ►IT! with 
total quark helicity —1/2. These two actually will contribute opposite sign to the matrix 
element in Eq. ()21|). because of the 75 in the operator (D L . Taking into account this fact, 
and substituting the above decomposition into Eq. (|21j). we find that the matrix element 
becomes 

i(ps z t\6 l \ps z1 )i - _ k {PS zi \6 L \PS- A )_ h 

-_i_(PS zt \6 L \PS zl )_i 2 + i_(PS z1 \6 l \PS z t)i . (24) 

It is easy to see that the above two lines are complex conjugates. In the following, we will 
consider the contribution from the first line, and the other one can be obtained immediately. 
For the subprocess TIT - *TIT) the contribution to the matrix element will be 

i_(PS^\6\PS xl ) h oc J ^(y[)^\y l ^(pl-ip\)T H (y i -y>- ViA _) , (25) 

where ip^ is the wave function for L z = Fock state and i/j^ for \L Z \ = 1 (their definitions 
are listed in the Appendix). Here we only show the contribution from the interference 
between ip^ and ip^ wave functions, and other interference contributions (e.g., the one 
with ifj^ and can be calculated similarly. Because proton is stable, the light-cone wave 
functions are real, i.e., (tp)* = ip. Meanwhile, for the IT! - ►IT! partonic process, we will 
have 

_i{PS xl \6\PS x i)_!_ * j i' {1 \y , ^ {3 \y i ,p l ±) ( - P? -ip\)T H {y i]y ' { ,p iL ) . (26) 

The hard partonic parts Tjj in the above two equations are identical to each other for the 
same diagram if we change all the quarks helicities. Thus we can sum their contributions 
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together, and the matrix element will be 

(PS x i\6\PS xi ) - (PS xl \6\PS x1 ) oc / ^ l \y[)^\y l1 p^){pl)T H {y l -y[-p li _) + h.c. (27) 



The linear expansion term of pi± from T# will be crucial to obtain nonzero contribution 
to the above matrix element when integrating over pij_. Otherwise, it will vanish. This 
expansion will introduce an additional suppression factor in (1 — x). For example, one of 
the propagators in Fig. 1 has the following expansion result, 

1 1 



(P3 - h) 2 (jfe-P ~h + P3±] 



Z/3^2± V 2/3 fc; 



2_L 



Substituting the above into Eq. (J27|) . and using the fact that p3j_ = — pu_ — p2±, we find 
that the above expansion will lead to a contribution as (1 — x) J d 2 pi±p^ ± (pi± ■ kj_) ip^ ~ 
(1 — x)A;^$ ( - 3 ' 4 ' ) , where <3>( 34 ) represents a combination of twist-four distribution amplitudes 
$4 and ^4, and the k\ factor will cancel out the same k\ in the denominator in Eq. (J21j) . 
This suppression feature applies to every propagator expansion containing the linear term 
of the intrinsic transverse momentum p±. Similarly, the Dirac wave function expansion in 
terms of p± will also be suppressed by (1 — x). 

The above analysis can be repeated for every diagrams in Fig. 2, and they contribute the 
same. So, the final power counting result for the TMD quark distribution Aq T will be, 

Ag T (x,fc ± )Ui ~ (1-x) 4 . (29) 

Similar analysis can be performed for the TMD quark distribution 5qL, which is defined 
through the following matrix element, 

Sq L (x,k ± ) = ¥* J ^^e-^- p+ +^ (PS z \W v (Ol + l^M\PS z ) . (30) 

If we choose Y = 7 X in the above equation, the TMD 6qi will become, 

SQL(x,k ± ) = ^J ^^e-^- p+ ^(PS z \O t \PS z ) , (31) 

where the operator Ot follows the definition in the subsection Sec. 11(b). In the above defini- 
tion, the proton is longitudinal polarized, and we can further write down explicitly in terms 
of the proton helicity states, 

6u(xM) = ^J ^^e"^ + +^^ . (32) 

Since the operator Ot is chiral-odd, it changes the quark helicity. However, in the above 
equation, we are calculating the hadron helicity conserved matrix elements, thus the nonzero 
quark orbital angular momentum projection must be taken into account in order to obtain 
nonzero results, as in the case of Ag^ i n the above. Following the above analysis, we find 
the power counting result for the TMD quark distribution Sqi, 

feU~(l-s) 4 , (33) 
which is one power of (1 — x) suppressed relative to the unpolarized quark distribution. 
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FIG. 3: Leading Feynman diagrams contributing to the naive time-reversal-odd TMD quark distri- 
butions at large-x. 



C. Naive time-reversal-odd quark distributions 

Now, we turn to the naive time-reversal-odd quark distributions. At leading order, we 
have two: the quark Sivers function and the Boer-Mulders function. The Sivers function 
represents the unpolarized quark distribution in a transversely polarized target, while the 
Boer-Mulders function represents the transversely polarized quark distribution in a unpolar- 
ized proton target. These two distributions are naive time- reversal-odd, and their existence 



require final state interactions 11]. The quark Sivers function is defined as, 



= J ^^^-^HPS\d\PS) , (34) 

where the operator O follows the above definition. Because the target is transversely polar- 
ized, again we will choose the x-direction for its polarization, and the Sivers function then 
becomes, 

Qt(x, k ± ) = J ^^e-^ 1 - ((PS*\6\PS Mi ) + (PS zl \6\PS z ,)) . (35) 

The above equation shows that the Sivers function is proportional to the matrix elements 
involving hadron helicity flip. Because the operator O is chiral-even, it conserves the quark 
helicities. To obtain the hadron helicity flip, we have to take into account the nucleon's 
light-cone wave function with nonzero quark orbital angular momentum, as in the previous 
two examples. 

Following the calculations in the last subsection for fc^-odd distribution Aqj-, we find 
that the Sivers function will depend on the following matrix element, 

l2 (PS z1 \d\PS zl ) l2 + _,(PS zl \6\PS^)^ 
+ UPS z i\6\PS zi ). + h {PS zi \6\PS z ^) h . (36) 
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Comparing with Eq. (j24jl . we find that only the sign changes in the above sum. This is 
because here we are probing the unpolarized quark, and we have to sum up different quark 
helicity contribution, while in Eq. ()24j) we are probing the longitudinal polarized quark and 
different quark helicity will contribute differently. From Eqs. (|25I26|) . we find that the final 
result for the above matrix element will be 

(PS x i\6\PS xL ) + (PS xl \d\PS^) ex J ^ 1) (y^ 3 \y l ) {-ip\)T H {y i]V ' { ,p iX ) +h.c. . (37) 

From this equation, we find that in order to generate a nonzero Sivers function, the hard 
scattering factor Th has to have an imaginary part. In a partonic hard scattering amplitude, 
the only imaginary part comes from the on- shell pole of some propagator. As we showed 
in the above analysis all the propagators in Fig. 2 are far off-shell except for the eikonal 
propagator from the gauge link. Thus, in order to obtain nonzero contribution to the Sivers 
function, we have to have eikonal propagator in the partonic Feynman diagrams. We have 
shown all leading order diagrams contribution to the Sivers function in Fig. 3, all of which 
have at least one eikonal propagator. For example, the eikonal propagator in Fig. 3(a) reads, 

77 1 r~ ■ = P— iir5{x - yi) , (38) 

n ■ [k — pi) + le x — y± 

where the first term is the principal value of the pole, and does not contribute to the Sivers 
function. Only the second term contribute to an imaginary part, which contains a delta 
function. This delta function will affect the power counting for the various factors in the 
evaluation of the matrix element of Eq. ()37)1 . This is because the delta function can be 
written as 5(x — y±) = 5 (y 2 + y$ — (1 — x)), which means that the variables y 2 and y 3 are 
limited to be order of (1 — x), i.e., y 2 ~ y% ~ 0(1 — x). All factors which depend on y 2 and y 3 
will have to be examined carefully to get the right power counting results. For example, in 
Fig. 3(a), the propagators at the left side of the cut will be affected by the above constraints. 
One of the gluon propagator reads, 

(39) 



y3j3(l- x ) 

because y% ~ (9(1— x) and (3 is order of unit. Unlike the case studied in the above subsections, 
this propagator does not lead to a suppression in (1 — x) for the Sivers function. Similarly, 
another gluon propagator and the quark propagator at the left side of the cut line are also 
finite at x — > 1. However, all the propagators at the right hand side of the cut line still scale 
as (1 — x), and the total four propagators there will contribute to a suppression factor of 
(1-x) 4 . 

Another consequence of this delta function is that the intrinsic p± expansion in the hard 
part has no additional suppression in (1 — x), which is very different from what we have 
in the last subsection for the fcj_-odd but naive time-reversal-even quark distributions. For 
example, in one of the above propagators, we can keep the intrinsic transverse momentum 
dependence and expand it up to the linear term, 

(40) 



(pj, - k 2 y -2k 2 ■ p 3 ~ -^k 2 2± - 2k 2L ■ p 3± (hi) V (kl) 
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The last equation of the above comes from the fact that 2/3 ~ o(l — x). After com- 
bining this expansion with the light-cone wave function, we get J p\ ± k 2 ± ■ p?,±^ 3 ^ oc 
k y ± y3& 4(2/1) or 2/2^4 (yi), where the k y ± factor will cancel that in the denominator in Eq. ()35|) . 

The final step in this analysis will be to factor out (1— x) from the wave function integral. 
Because the variables y 2 and 2/3 are constrained to be order of (1 — x) in the x — > 1 limit, we 
have to examine the wave function integral where the end-point behavior of the light-cone 
wave function will be important. For example, one contribution in the above analysis is the 
integral of 2/3^4(2/1, 2/2; 2/3)- We can factor out the overall dependence on (1 — x) from the 
integral, 

dy 1 dy 2 dy 3 5(l - y x - y 2 - 2/3) %i - x)y 3 ^ 4 (yi, 2/2, Vs) 
= (1 - x) 3 J d(dC5(l - C - COCCI - x)**(x, C(l - x), C'(l - x))/{\ - xf , (41) 

where we have re-parameterized y 2 = C(l — x) and 2/3 = £'(1 — x )- The above integral 
depends on the end-point behavior of the twist-four distribution amplitudes. As we showed 
in the Appendix, the twist-four distribution amplitude behaviors as 2/3^4(2/1, 2/2, 2/3) °^ 2/12/22/3 
at the end-point region. From this, we find that 

C(l-x)$ 4 (x,£(l- x), C'(l - x)) 
x-^i (1 - x) 



lim ^^-^K^-^^--JJ = finite (42) 



Thus the wave function integral indeed contains a suppression factor (1 — x) 3 . In addition, 
the Dirac matrix traces will also result into a power dependence on (1 — x). This can 
be calculated straightforwardly, and we find the Dirac traces from Fig. 3(a) contribute to 
a power term as (1 — x)~ 2 . By summarizing the power counting results from the above 
analysis and also taking into account the phase spaces integral factor (1 — x) _1 in Eq. 
we find the Sivers function will have the following power behavior, 

Qt(x, fcjJU-* oc (1 -x) 4 , (43) 

which is (1 — x) suppressed relative to the unpolarized quark distribution. Similar calcula- 
tions can be performed for all other diagrams in Fig. 3, and they all contribute to a power 
behavior of (1 — x) 4 for the quark Sivers function. 

The same analysis can be performed for another naive time-reversal odd distribution, the 
so-called Boer-Mulders function Sq, which is defined as 

6q = I Wf 6 '^^^ < P l^(OT + 7 i 7 5 ^(0)|P> . (44) 

If we choose 7* = 7 X in the above equation, the TMD 5q will become, 

8q(x, k±) = ^J ^e~^- p+ ^(P\O t \P) , (45) 

where the operator O t follows the definition in the above. Unlike the TMD quark distribution 
SqL, in the above definition the proton is unpolarized. Thus the explicit expression for 5q in 
the proton helicity states is 

8q(x,k ± ) = ^J ^y^- txi ' P++l ^\ ((PS z ,\6 t \PS z ,) + (PS zl \6 t \PS zl )) . (46) 
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Comparing this expression to Eq. in the last subsection, we find that 5q depends on 
the sum of the two matrix elements, whereas for 8qi it is the difference. As in the analysis 
for the Sivers function, we find that we need an imaginary part from the hard part, and the 
same set of diagrams in Fig. 3 contribute. The final result for its power counting will be 

Sq\ x ^ ~ (1 - x) 4 , (47) 

which is again one power of (1 — x) suppressed relative to the unpolarized quark distribution. 

D. Comparison with the Power Counting for the GPD E 

Summarizing the results in the last two subsections, we find that the &j_-odd TMD 
quark distributions are suppressed by a relative factor — to the fcj_-even ones (e.g., 
the unpolarized quark distributions). As we mentioned in the introduction, this can also be 
understood by the interpretation of these TMD quark distributions in terms of the overlaps 
of the light-cone wave functions of L z = and \L Z \ = 1 Fock states. 

As in the case of the Pauli form factor, the generalized parton distribution (GPD) E 
and the Sivers function gy involve the overlap of initial- and final-state light-front-wave- 
functions (LFWFS) which differ by one unit of orbital angular momentum |37[. In contrast 
to the Sivers function which is suppressed by one power of (1 — x) relative to the unpolarized 
distribution, one finds the GPD E falls as two-powers (1 — x) 2 faster than the spin-conserving 
GPD H at large-x [38]. The power of (1 — x) n thus differs when we compare the E GPD 
arising in spin-flip deeply virtual Compton scattering (DVCS) -y*Pi ~* IPi an d the Sivers 
function qx arising in polarized electroproduction 7*pj — > ttX. In the following, we will 
briefly comment why this happens. 

It is useful to use the symmetric light-front (LF) frame where the transverse momenta 
of the initial and final state proton momentum changes from p™ ltial = (p ± — |Aj_) to 
p final _ _|_ |A_l). The struck quark in DVCS is evaluated at k± + ~(1 — x)A^ in the 

final-state LFWF and k± — |(1 — x)Aj_ in the initial-state LFWF, as in the Drell-Yan-West 
(DYW) formula for current matrix elements (3^. 

The E GPD requires evaluating the spin-flip deeply virtual Compton amplitude which 
is linear in the transverse momentum transfer to the proton Aj_. This kinematic factor 
arises from the extra angular momentum of the initial- or final-state LFWF with argument 
±|(1— x)A±. In addition, the orbital angular momentum dynamics of the LFWF introduces 
a factor of (1 — x). Thus E ~ (1 — x) 2 H as x — > 1. 

In contrast, when we evaluate the Sivers SSA for SIDIS 7*p — > up', the dynamics of the 
orbital angular momentum in the LFWF gets expressed as the transverse momentum p n ± 
of the produced pion, not the change in the transverse momentum A^ of the proton. Thus 
the second factor of (1 — x) does not appear in the Sivers function. We thus have the power 
counting rule: E ~ (1 — x)qr ~ (1 — x) 2 H as x — > 1. 

E. Power counting for the integrated quark distributions at leading and higher- 
twist 

From the power counting results for the TMD quark distributions in the last subsections, 
we can further derive the power counting rule for the integrated quark distributions at large- 
x when integrating over the transverse momentum. For example, the integrated unpolarized 
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quark distribution can be written as 



q(x) = j d 2 k ± q(x, k±) . (48) 

Similar equations also hold for the longitudinal polarized quark distribution and transversity 
quark distribution. From the power counting of these relevant TMD quark distributions, 
we can immediately see that the integrated quark distributions have the following power 
counting rule at x — > 1, 

q(x) ~ (1 -a;) 3 , Aq L (x) ~ (1 - x) 3 , 5q T (x) ~ (1 - xf . (49) 

To obtain the above power counting results for the integrated quark distributions, we have 
assumed that the /c^-integral decouples from the x— distributions of the partons [i^, 
Although the upper limit of the k± integral might depend on (1 — x), the bulk of this 
integration comes from the lower bound, which will not affect the (1 —x) power counting for 
the integrated parton distributions [29j. The latest comparison of the above power counting 
predictions with experiment can be found in |4p| . 

The /c_i_-moment of the fcj_-odd TMD quark distributions are related to the twist-three 
parton distributions. For example, the twist-three parton distribution gr(x) is related to 
the TMD quark distribution Ag T [3Cj, 



9t(x) = j d 2 k ± kiAq T (x, k ± ) , (50) 



h L {x) = — / d 2 k ± ki5q L {x, k±) • (51) 



and for h^x), 



Of course, caution has to be taken when we apply the above equations j30|. From the power 
counting rule for the relevant TMD quark distributions, we find the following power behavior 
for these two twist-three parton distributions, 

g T (x)~(l-x)\ h L {x)~{l-x) 4 . (52) 

/c_i_-moment of the naive-time-reversal-odd TMD quark distributions are also related to the 
twist-three parton distributions, which have been shown in literature wJl, for example, 



T F (x) = — J d k ± k ± q T (x,k ± ) , (53) 

where Tp is the so-called Qiu-Sterman matrix element , and is responsible to the SSA for 
inclusive hadron production in hadronic collisions. Similarly, the /c^-moment of Sq corre- 
sponds, 

where Tp is defined as 



TP(x) = — j d 2 k ± ki5q(x } k ± ) , (54) 



T { ;\x u x 2 ) = J d <^L e ^P + v- e i^-^C- (p\ ma +"g F + a (C)Hv-)\P) , (55) 
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(x, k±) 




(e) (4 

FIG. 4: Feynman diagrams contribution to the transverse-momentum-dependent quark distribution 
in Pion at large-x. 

and Tp\x) = Tp(x,x). From the power counting rule of the relevant naive-time- reversal- 
odd TMD quark distributions, we find the power counting rule for these two twist-three 
parton distributions, 

T F (x) ~ (1 - x)\ T { p\x) ~ (1-x) 4 , (56) 
which are (1 — x) suppressed relative to the unpolarized quark distribution. 



F. TMD quark distributions in Pion 

A similar analysis can be carried out for the TMD quark distributions of pion. Because 
the pion is a spin-0 particle, there are only two leading-order TMD quark distributions: 
the unpolarized quark distribution q 7r (x,k±) and the Boer-Mulders function Sq 7r (x,k±). At 
large-x, their dependence on x can be calculated from the diagrams shown in Fig. 4. The un- 
polarized quark distribution at large-x will have contributions from all of the four diagrams, 

u«(x, k ± ) = -4^(1 - xfo? s C F / ——^{zx)^(ztfT H (zx, z 2 ) , (57) 

J Zl z 2 

where f w is the decay constant of pion, and <& % {z) is the leading-twist quark distribution 
amplitude. Obviously, the quark distribution has (1 — x) 2 power behavior at large-x. This 
is consistent with the Gribov-Liptov relation . 

For the Boer-Mulders function, because it is naive-time-reversal-odd, we have to take into 
account the interference between L z = and \L Z \ = 1 Fock states of pion light-cone wave 
functions, and also the gauge link is important to obtain a phase difference. Following the 
same analysis in the previous subsections, we find that the Boer-Mulders function of pion 
has the same power counting result as the unpolarized quark distribution, 

Sq«(x,k ± ) ~ (1-x) 2 . (58) 

These two distribution functions having the same power behavior at large-x, is not because 
the Boer-Mulders function of pion gets enhancement, but because the unpolarized quark 
distribution of pion is suppressed by one power of fl — x) compared to the usual power 
counting results for parton distributions of hadrons 24[ 2?] |. 
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III. CONCLUSION 



In this paper, we have performed a perturbative analysis of the transverse momentum 
dependent quark distributions at large x. A generalized power counting rule has been de- 
rived for the leading order TMD quark distributions, and we have found that the k±-even 
distributions all scale as (1 — x) 3 , whereas the fcj_-odd ones as (1 — x) 4 , including the naive 
time-reversal-even and -odd distributions. In particular, we have shown that the quark 
Sivers function has power behavior of (1 — x) 4 , which is (1 — x) suppressed relative to the 
unpolarized quark distribution. For the TMD quark distributions of pion, we find that the 
Boer-Mulders function has the same power behavior as the unpolarized quark distribution, 
scaling as (1 — x) 2 in the limit. These results provide important guidelines for the parame- 
terizations of the transverse momentum dependent parton distributions and the quark-gluon 
correlation functions in the phenomenological studies. 

In our analysis we have not included the effects of perturbative QCD evolution. In 
fact, the evolution of parton distributions at large x with photon virtuality Q 2 is suppressed 
compared to the usual DGLAP evolution |25( because the struck quark is a bound constituent 
of the target hadron. In particular, in the limit of x — > 1, the virtuality of the struck quark 
becomes highly space-like, and evolution is effectively quenched. Thus the power counting 
of structure functions at large x is not affected by evolution |25| , allowing duality with the 
power-law falloff of exclusive channels at fixed W 2 . Another important point has to be kept 
in mind is the large logarithms associated with the parton distributions in the x — > 1 limit, 
in terms of a™log m (l/(l - x)) for m < 2n 0, HI, Hi, 0. All these effects will of course 
introduce additional theoretical uncertainties when we apply the power counting rule to the 
parton distributions at large-x. 

Acknowledgments 

We thank Xiangdong Ji, Jianwei Qiu, and Werner Vogelsang for their comments. S.J.B. 
is supported by the Department of Energy under contract number DE-AC02-76SF00515, 
F. Y. is grateful to RIKEN, Brookhaven National Laboratory and the U.S. Department 
of Energy (contract number DE-AC02-98CH10886) for providing the facilities essential for 
the completion of their work. Finally, F.Y. thanks Beijing Institute of Modern Physics, 
Institute of Theoretical Physics, Beijing University, for its warm host when part of this 
work was finished. 



17 



APPENDIX A: LIGHT-CONE WAVE FUNCTIONS AND DISTRIBUTION AM- 
PLITUDES OF NUCLEON 



In this Appendix we list the light-cone wave functions for the three-quark Fock states of 
nucleon as reference |30| . 



\ps z1 ) 



d[l]d[2]d[3] |^ (1) (l,2,3)u{(l) (u\ (2)d{(3) - d\(2)u\(3)) |0> 

+ ((p£ + ipl)ft 3) (h 2, 3) + (p x 2 + ipl)^\l, 2, 3)) 
(«{(1)4(2)4(3) -4(1)4(2)^(3)) |0)} , 

J d[l]d[2]d[3] 2, 3)4(1) (4(2)4(3) - 4(2)4(3)) |0) 

+ ((P? -^ (3) (1,2,3) + (p§ -^)^ (4) (1,2,3)) 
4(1)4( 2 )4( 3 ) — rf{(l)^{(2)w{(3)) |0)| , (Al) 



where the argument i is the shorthand for quark momentum variables yi and pi±, and the 
measure for the quark momentum integrations is 



d[l]d[2]d[3] = V2 



dy 1 dy 2 dy 3 d p 1± d p 2± d p 3± 



V2yi2y 2 2y 3 (2vr)9 
x2tt<5(1 -yi - y 2 - y 3 )(2ir) 2 5 {2) (p 1± + p 2± + p 3± ) . 



(A2) 



^(1,3,4) 

are the light-cone wave function amplitudes for the three quark Fock state expansion 
of nucleon. tp^ corresponds to the L z = Fock state component, and ip^ 3 ' 4 ^ for \L Z \ = 1 ones. 
These light-cone wave functions were used in our analysis for the large-x quark distributions. 

In order to get Eq. (131, th e light-cone wave functions have to be converted into the quark 
distribution amplitudes 25||. For example, we can integrate out the transverse momentum 
in the leading Fock state light-cone wave function, and define the twist-three amplitude, 



s s (w) =2V6 / d ^ l±d ^f^ ^ (p 1± + m + ns)^ (1,2,3) 



(A3) 



For \L Z \ = 1 states, we have to keep linear term in the p± expansion of the hard factor, and 
combine them with the light-cone wave function, which will lead to the twist-four distribution 
amplitudes of the nucleon 

0,0, 



^4(2/1,2/2,2/3) 



$4(2/2,2/1,2/3) 



2a/6 f d 2 p xl _d 2 p 21 _d 2 p 3 ^ 



2/2M 

XP*2_L 

2^6 
2/3~M 

XP3± 



5 {2) (P1± +P2±+P3±) 



(2tt) 6 

p 1± ^ (3) (l,2,3)+p 2± ^ 4 )(l,2,3) 

d 2 Pi±d 2 p 2± d 2 p 3± (2) 

^ S { , (Pi±+P2±+P3l) 

p 1± ^ 3 \l,2,3)+p 2± ^ 4 \l,2,3) 



(A4) 
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The explicit expressions for these distribution amplitudes are not necessary for the power- 
counting analysis. However, the end-point behavior at j/j — > 1 is needed for the power 
counting of the naive time-reversal-odd TMD quark distributions. We note that in the end- 
point region, these distribution amplitudes have the following behaviors: $3(2/4) oc 2/12/22/3, 
2/2^/4(2/4) oc 2/12/22/3? an d 2/3^4 (Vi) oc 2/12/22/3 HI- From this, we immediately find that the 



end-point behavior of the p^-moment of the light-cone wave functions. For example, 

J rfVl±^V2±^V3±5 (2) (Pl-L + P2± + P3±){Pi± ■ Pj±W 3A \yi,y2, 2/3) lend point ~ 2/l2/22/3 , (A5) 

where i, j = 1, 2, 3. These properties have been used in our analysis. 



[1] see for example: G. Bunce et al., Phys. Rev. Lett. 36, 1113 (1976); D. L. Adams et al. [E581 
and E704 Collaborations], Phys. Lett. B 261, 201 (1991); D. L. Adams et al. [FNAL-E704 
Collaboration], Phys. Lett. B 264, 462 (1991); K. Krueger et al, Phys. Lett. B 459, 412 
(1999). 

[2] A. Airapetian et al. [HERMES Collaboration], Phys. Rev. Lett. 84, 4047 (2000); A. Airapetian 
et al. [HERMES Collaboration], Phys. Rev. D 64, 097101 (2001); Phys. Rev. Lett. 94, 012002 
(2005); M. Diefenthaler [HERMES Collaboration], arXiv:hep-ex/0507013| 

[3] V. Y. Alexakhin et al. [COMPASS Collaboration], Phys. Rev. Lett. 94, 202002 (2005). 

[4] H. Avakian et al. [CLAS Collaboration], Phys. Rev. D 69, 112004 (2004) 
|arXiv:hep-ex70 301005 ; H. Avakian, P. Bosted, V. Burkert and L. Elouadrhiri [CLAS 
Collaboration], AIP Conf. Proc. 792, 945 (2005) |arXiv:nucl-ex/0509032| . 

[5] J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 92, 171801 (2004) 
|arXiv:hep-ex/0310058| . 

[6] S. S. Adler [PHENIX Collaboration], Phys. Rev. Lett. 95, 202001 (2005) 
|arXiv:hep-ex/0507073| . 

[7] F. Videbaek [BRAH MS Collaboration], AIP Conf. Proc. 792, 993 (2005); 
|arXiv:nucl- ex/0601008 J. H. Lee [BRAHMS Collaboration], talk presented at the "14th 
International Workshop on Deep Inelastic Scattering (DIS 2006)", Tsukuba, Japan, April 
20-24, 2006. 

[8] A. V. Efremov and O. V. Teryaev, Sov. J. Nucl. Phys. 36, 140 (1982) [Yad. Fiz. 36, 242 

(1982) ]; A. V. Efremov and O. V. Teryaev, Phys. Lett. B 150, 383 (1985). 

[9] JW. Qiu and G. Sterman, Phys. Rev. Lett. 67, 2264 (1991); Nucl. Phys. B 378, 52 (1992); 

Phys. Rev. D 59, 014004 (1999). 
[10] D. W. Sivers, Phys. Rev. D 41, 83 (1990); Phys. Rev. D 43, 261 (1991). 
[11] S. J. Brodsky, D. S. Hwang and I. Schmidt, Phys. Lett. B 530, 99 (2002); Nucl. Phys. B 642, 

344 (2002). 

[12] J. C. Collins, Phys. Lett. B 536, 43 (2002). 

[13] X. Ji and F. Yuan, Phys. Lett. B 543, 66 (2002); A. V. Belitsky, X. Ji and F. Yuan, Nucl. 

Phys. B 656, 165 (2003). 
[14] D. Boer, P. J. Mulders and F. Pijlman, Nucl. Phys. B 667, 201 (2003). 

[15] J. C. Collins and D. E. Soper, Nucl. Phys. B 193, 381 (1981) [Erratum-ibid. B 213, 545 

(1983) ]; Nucl. Phys. B 197, 446 (1982); J. C. Collins and D. E. Soper, Nucl. Phys. B 194, 
445 (1982). 



19 



[16] J. C. Collins, D. E. Soper and G. Sterman, Nucl. Phys. B 250, 199 (1985). 

[17] X. Ji, J. P. Ma and F. Yuan, Phys. Rev. D 71, 034005 (2005); Phys. Lett. B 597, 299 (2004); 

JHEP 0507, 020 (2005) |arXiv:hep-ph/0503015| . 
[18] J. C. Collins and A. Metz, Phys. Rev. Lett. 93, 252001 (2004). 

[19] X. Ji, J. W. Qiu, W. Vogelsang and F. Yuan, Phys. Rev. Lett. 97, 082002 (2006) 



arXiv:hep-ph/0602239 ; Phys. Rev. D 73, 094017 (2006) arXiv:hep-ph/0604023 ; Phys. Lett. 
B 638, 178 (2006) |arXiv:hep-ph/0604128| . 
[20] M. Anselmino, M. Boglione, U. D'Alesio, A. Kotzinian, F. Murgia and A. Prokudin, Phys. 
Rev. D 72, 094007 (2005) [Erratum-ibid. D 72, 099903 (2005)] |arXiv:hep-ph/0507181| ; Phys. 
Rev. D 71, 074006 (2005) |arXiv:hep-ph/0501196| . 
[21] W. Vogelsang and F. Yuan, Phys. Rev. D 72, 054028 (2005) |arXiv:hep-ph/0507266| . 
[22] J. C. Collins, A. V. Efremov, K. Goeke, S. Menzel, A. Metz and P. Schweitzer, Phys. Rev. 

J. C. Collins et al, Phys. Rev. D 73, 094023 



014021 (2006) 


arXiv:hep-ph/0509076 


arXiv:hep-ph/0511272 . 



[23] C. Kouvaris, J. W. Qiu, W. Vogelsang and F. Yuan, arXiv:hep-ph/0609238 
[24] G. R. Farrar and D. R. Jackson, Phys. Rev. Lett. 35, 1416 (1975). 
[25] G. P. Lepage and S. J. Brodsky, Phys. Rev. D 22, 2157 (1980). 
[26] A. H. Mueller, Phys. Rept. 73, 237 (1981). 

[27] J. F. Gunion, P. Nason and R. Blankenbecler, Phys. Rev. D 29, 2491 (1984). 

[28] S. J. Brodsky, M. Burkardt and I. Schmidt, Nucl. Phys. B 441, 197 (1995) 

|arXiv:hep-ph/9401328] . 
[29] X. Ji, J. P. Ma and F. Yuan, Phys. Lett. B 610, 247 (2005) |arXiv:hep-ph/0411382| . 
[30] X. Ji, J. P. Ma and F. Yuan, Nucl. Phys. B 652, 383 (2003) |arXiv:hep-ph/0210430| . 
[31] S. J. Brodsky and S. Gardner, arXiv:hep-ph /0608219| 
[32] S. J. Brodsky and S. D. Drell, Phys. Rev. D 22, 2236 (1980). 

[33] P. J. Mulders and R. D. Tangerman, Nucl. Phys. B 461, 197 (1996) [Erratum-ibid. B 484, 
538 (1997)] larXiv:hep-ph/9510301| ; D. Boer and P. J. Mulders, Phys. Rev. D 57, 5780 (1998) 
|arXiv:hep-ph/9711485| . 

[34] R. L. Jaffe and X. Ji, Nucl. Phys. B 375, 527 (1992). 

[35] J. C. Collins, Nucl. Phys. B 396, 161 (1993) |arXiv:hep-p h/9208213] . 

[36] J. Soffer, Phys. Rev. Lett. 74, 1292 (1995) |arXiv:hep-ph/9409254| . 

[37] S. J. Brodsky, M. Diehl and D. S. Hwang, Nucl. Phys. B 596, 99 (2001) 

|arXiv:hep-ph /0009254] . 
[38] F. Yuan, Phys. Rev. D 69, 051501 (2004) |arXiv:hep-ph/0311288| . 

[39] S. D. Drell and T. M. Yan, Phys. Rev. Lett. 24, 855 (1970); G. B. West, Phys. Rev. Lett. 24, 
1206 (1970). 

[40] X. Zheng et al. [Jefferson Lab Hall A Collaboration], Phys. Rev. Lett. 92, 012004 (2004) 

|arXiv:nucl-ex/0308011| ; Phys. Rev. C 70, 065207 (2004) |arXiv:nucl-ex/0405006| . 
[41] V. N. Gribov and L. N. Lipatov, Sov. J. Nucl. Phys. 15, 675 (1972) [Yad. Fiz. 15, 1218 (1972)]. 
[42] G. P. Korchemsky, Mod. Phys. Lett. A 4, 1257 (1989). 
[43] G. Sterman, Nucl. Phys. B 281, 310 (1987). 

[44] V. Braun, R. J. Fries, N. Mahnke and E. Stein, Nucl. Phys. B 589, 381 (2000) [Erratum-ibid. 

B 607, 433 (2001)] |arXiv:hep-ph/0007279| . 
[45] A. V. Belitsky, X. Ji and F. Yuan, Phys. Rev. Lett. 91, 092003 (2003) |arXiv:hep-ph/0212351| . 



20 



